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We investigate the dynamical one-armed spiral instability in differentially rotating stars with both eigcnmodc 
analysis and hydrodynamic simulations in Newtonian gravity. We find that the one-armed spiral mode is 
generated around the corotation radius of the star, and the distribution of angular momentum shifts inwards 
the corotation radius during the growth of one-armed spiral mode. We also find by investigating the distribution 
of the canonical angular momentum density that the low T/[W] dynamical instability for both m = 1 and m = 2 
mode, where T is the rotational kinetic energy and W is the gravitational potential energy, is generated around 
the corotation point. Finally, we discuss the feature of gravitational waves generated from these modes. 



1. Introduction 

Stars in nature are usually rotating and may be sub- 
ject to non-axisymmetric rotational instabilities. An 
exact treatment of these instabilities exists only for in- 
compressible equilibrium fluids in Newtonian gravity, 
e -g- ..3,0- For these configurations, global rotational 
instabilities may arise from non-radial toroidal modes 
e lm{ p (where rn = ±1, ±2, . . . and ip is the azimuthal 
angle) . 

For sufficiently rapid rotation, the m = 2 bar mode 
becomes either secularly or dynamically unstable. The 
onset of instability can typically be marked by a crit- 
ical value of the dimensionless parameter (3 = T/|W|, 
where T is the rotational kinetic energy and W the 
gravitational potential energy Uniformly rotating, in- 
compressible stars in Newtonian theory are secularly 
unstable to bar-mode formation when fj > /3 SGC ~ 
0.14. This instability can grow only in the presence of 
some dissipative mechanism, like viscosity or gravita- 
tional radiation, and the associated growth timescale 
is the dissipative timescale, which is usually much 
longer than the dynamical timescale of the system. 
By contrast, a dynamical instability to bar-mode for- 
mation sets in when (3 > /?dyn — 0.27. This instability 
is independent of any dissipative mechanism, and the 
growth time is the hydrodynamic timescale. 

There are several papers indicating that dynamical 
instability of the rotating stars occurs at low T/|W|, 
which is far below the standard criterion of dynamical 
instability in Newtonian gravity. Tohline and Hachisu 
Q find in the self-gravitating ring and tori that m = 2 
dynamical instability occurs around T/|W| ~ 0.16 in 
the lowest case in the centrally condensed configu- 
rations. Shibata, Karino, and Eriguchi |j| find that 
m = 2 dynamical instability occurs around T/|W| ~ 
O(10" 2 ) in the high degree (n c /n eq rj 10) of differen- 
tial rotation. Note that f2 c and f2 oq are the angular 
velocity at the center and at the equatorial surface, re- 
spectively. Centrella et al. 5] found dynamical m = 1 
instability around T/|W| ~ 0.09 in the n = 3.33 poly- 
tropic toroidal star with high degree (f2 c /f2 cq = 26) 



of differential rotation, and Saijo, Baumgarte, and 
Shapiro [|J extended the results of dynamical m = 1 
instability to n > 2.5, fi c /f2 e q ^ 10. Note that n is 
the polytropic index of the star. 

There are some indications that corotation reso- 
nance triggers dynamical bar instability. Papaloizou 
and Pringle found that non-selfgravitating tori 
with constant specific angular momentum are unsta- 
ble to low order non-axisymmetric modes and that 
the modes grow on a dynamical time-scale. Watts, 
Andersson, and Jones [g argue that the shear insta- 
bilities occur when the degree of differential rotation 
exceeds a critical value and the /-mode develops a 
corotation point associated with the presence of a con- 
tinuous spectrum. They also point out that dynamical 
bar instability that Shibata et al. found is in the 
corotation band. 

Our purpose in this paper is to investigate the na- 
ture of dynamical m = 1 instability with both eigen- 
modc analysis and hydrodynamical analysis. A non- 
linear hydrodynamical simulation is indispensable for 
investigation of evolutionary process and final out- 
come of instability. The nature of instability as a 
source of gravitational wave, which interests us most, 
is only accessible through non-linear hydrodynamical 
computations. On the other hand, a linear eigenmode 
analysis is in general easier to approach the dynamical 
instability of given equilibria and it may be helpful to 
have physical insight on the mechanism and the ori- 
gin of instability. Therefore a linear eigenmode analy- 
sis and a non-linear simulation are complementary to 
each other and they both help us to understand the 
nature of dynamical instability. 

For a hydrodynamical simulation, we used the nu- 
merical code developed in Ref. |6|, while we intro- 
duced a toy cylinder model that mimics differentially 
rotating stars to study the instability. Self-gravitating 
cylinder models have been used to study general dy- 
namical nature of such gaseous masses as stars, accre- 
tion disks and of stellar system as galaxies. Although 
there is no cylinder with infinite length in reality, it 
is expected to share some qualitative similarities with 
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Table I Three differentially rotating equilibrium stars that trigger dynamical instability 



Model 


n 


d/ Rcq b 


Rp / Rcq 


^c/^cq 


pc 1 pmax 


^maxd / -^cq 


f T/\W\ 3 


dominant unstable mode 


I 


3.33 


0.20 


0.413 


26.0 


0.491 


0.198 


0.146 


m — 1 


II 


1.00 


0.20 


0.250 


26.0 


0.160 


0.383 


0.119 


m = 2 


III 


1.00 


1.00 


0.250 


2.0 


0.837 


0.388 


0.277 


m = 2 



a n: Polytropic index 
'-Req: Equatorial radius 
c Rpl'. Polar radius 

d Q c : Central angular velocity; f2 eq : Equatorial angular velocity 
e p c : Central density; p max : Maximum density 
f R max( j '■ Radius of maximum density 

9 T: Rotational kinetic energy; W: Gravitational potential energy 



realistic astrophysical objects, e.g. Q. In fact, it has 
served as a useful model to investigate secular and 
dynamical instabilities of rotating masses. 

This paper is organized as follows. In §[2]we present 
our hydrodynamical results of dynamical one-armed 
spiral and dynamical bar instability. We present our 
diagnosis of dynamical m = 1 and m = 2 instability 
by using canonical angular momentum in § and 
briefly summarize our findings in § 0] Throughout 
this paper we use gravitational units with G = 1. A 
more detailed discussion will be presented in Ref. . 



2. Dynamical instabilities in differentially 
rotating stars 

First we explain features of our initial data sets of 
differentially rotating stars on which we performed 
non-linear hydrodynamical computations. We assume 
a polytropic equation of state, 



P 



up 



l + l/r. 



(1) 



where P is the pressure, K is a constant, p is the den- 
sity, n is the polytropic index. One feature of the poly- 
tropic equation of state is that all matter quantities 
can be renormalized in terms of k so that n does not 
explicitly appear. We also assume the "j-constant" 
rotation law as 



Jo 



d 2 



where fi is the angular velocity, jo is a constant pa- 
rameter with units of specific angular momentum, and 
w is the cylindrical radius. The parameter d deter- 
mines the length scale over which SI changes; uniform 
rotation is achieved in the limit d — > oo, with keep- 
ing the ratio_?'o/d 2 finite. We choose the same data 
sets as Ref. Q for investigating low T/|W| dynamical 
instability in differentially rotating stars. We also con- 
struct the equilibrium of a star with weak differential 
rotation in high T/|W|, which excites the standard 



dynamical bar instability, e.g. Q. The characteristic 
parameters are summarized in Table [I] 

To enhance any m = 1 or m = 2 instability, we 
disturb the initial equilibrium density p cq by a non- 
axisymmetric perturbation according to 



\ -Kcq 



R 2 



(3) 



where R cq is the equatorial radius, with <jW = « 
1.7 — 2.8 x 10~ 3 in all our simulations. We also intro- 
duce "dipole" D and "quadrupole" Q diagnostics to 
monitor the development of m — 1 and m = 2 modes 
as 



D = (e mllp ) , = 

\ / rn—1 



l 

M 



x + iy 3 



\A 2 + V' 



(fx, 



(4) 



« - \ e ) m =2 ~ M J P x 2 +y 2 ' 

(5) 

respectively. 

We compute approximate gravitational waveforms 
by using the quadrupole formula. In the radia- 
tion zone, gravitational waves can be described by 
a transverse-traceless, perturbed metric hf^ with re- 
spect to a flat spacetime. In the quadrupole formula, 
h^ T is found from 



■TT 



2 d 

rM 2 13 ' 



(G) 



(2) where r is the distance to the source, ly the 
quadrupole moment of the mass distribution, and 
where TT denotes the transverse-traceless projection. 
Choosing the direction of the wave propagation to be 
along the rotational axis (z-axis), the two polarization 
modes of gravitational waves can be determined from 



h+^\(h T T -h T y T) and h x ^h^. (7) 

For observers along the rotation axis, we thus have 
rh+ 1 d 2 



TT 



M 2M dt 2 ^ xx 



I TT ) 



(8) 
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Figure 1: Diagnostics D and Q as a function of t/P c for 
3 differentially rotating stars (see Table HJ. Solid and 
dotted lines denote the values of D and Q, respectively. 
The Roman numeral in each panel corresponds to the 
model of the differentially rotating stars, respectively. 
Hereafter P c represents the central rotation period. 



M 



1 d -fTT 

Mdt^ xv ' 



(9) 



The time evolutions of dipole diagnostic and the 
quadrupole diagnostic are plotted in Fig. We de- 
termine that the system is stable to m = 1 (m = 2) 
mode when the dipole (quadrupole) diagnostic re- 
mains small throughout the evolution, while the sys- 
tem is unstable when the diagnostic grows exponen- 
tially at the early stage of the evolution. It is clearly 
seen in Fig. ^that the star is more unstable to the 
one-armed spiral mode for model I, and more unstable 
to the bar mode for models II and III. In fact, both 
diagnostics grow for model I. The dipole diagnostic, 
however, grows more markedly than the quadrupole 
diagnostic, indicating that the m = 1 mode is the 
dominant unstable mode. 
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Figure 2: Density contours in the equatorial plane for 3 
differentially rotating stars (see Table EJ. Models I, II, 
and III are plotted at the parameter (t/P c , Pmax/ PmLi) = 
(16.3, 3.63), (132, 1.25), and (5.49, 1.20), where p max is 
the maximum density, Pmax is the maximum density at 
t = 0, and R is the equatorial radius at t = 0. The 
contour lines denote densities 

p/pmax = 10- (1(i ~ i)d (i = 1, ■ • • , 15) for model I and 
p/p m ax = 6.67(16 - i) X 10-' 2 (i = 1, • •• , 15) for models II 
and III, respectively. 
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The density contour of the differentially rotating 
stars are shown in Fig. [3 ft is clearly seen in Fig. [21 
that one-armed spiral structure is formed at the inter- 
mediate stage of the evolution for model I, and that 
bar structure is formed for models II and III once the 
dynamical instability sets in. 

We also show gravitational waves generated from 
dynamical one-armed spiral and bar instability in Fig. 
13 For to = 1 modes, the gravitational radiation 
is emitted not by the primary mode itself, but by 
the m = 2 secondary harmonic which is simultane- 
ously excited, albeit at a lower amplitude. Unlike the 
case for bar-unstable stars, the gravitational wave sig- 
nal does not persist for many periods, but instead is 
damped fairly rapidly. 



3. Canonical Angular Momentum 

We introduce canonical angular momentum follow- 
ing Ref. 01 t° diagnose oscillations in rotating fluids. 
For adiabatic linear perturbations on a perfect fluid 
configuration in stationary, axisymmetric spacetime, 
it is possible to introduce canonical conserved quanti- 
ties. Since we only use canonical angular momentum 
J c in this paper, we write down its explicit form as 

J c = 777 J (3?[cr] - mfl)p\Z\ 2 d 3 x 

-2m J pmfl ■ ^[^^*]d 3 x, (10) 

where a is the eigenfrequency, £ l is Lagrangian dis- 
placement vector. Note that total canonical angular 
momentum becomes zero when dynamical instability 
sets in. 

Next we apply the method of canonical angular mo- 
mentum to the linearized oscillations of a cylinder. We 
prepare two m = 1 stable modes (A-i, A-iii) and one 
to = 1 unstable mode (A-ii), summarized in Table ITTI 
We plot the integrand of canonical angular momen- 
tum wjc 

j c = 777(5R[cr] - mQ)p\£\ 2 - 2mpzun ■ (11) 

for 777 = 1 mode in Fig. 0] We define corotation radius 
H7 cr t of modes as 5ft[cr] — mCl(zo) = 0. This means that 
the pattern speed of mode coincides with the local ro- 
tational frequency of background flow there. We find 
that the distribution of canonical angular momentum 
changes its sign around corotation radius in m = 1 
unstable case, and that it is positive for no < m CI t, 
while it is negative for w > vj CIt - The behavior of the 
canonical angular momentum suggests us that this in- 
stability is related to the existence of corotation point 
inside the cylinder. We also find that the positive 
case of the canonical angular momentum (A-iii) cor- 
responds to the case when the pattern speed of the 
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Figure 3: Gravitational waveform for 3 differentially 
rotating stars (see Table as seen by a distant observer 
located on the rotational axis of the equilibrium star. 

Table II Parameters for equilibrium gaseous fluid and 
eigenfrequency 

Model flc/q h T/\W\ a/Q c 1 CTcrt/gg 3 

A-i 11.34 0.460 -0.245 

A-ii 11.34 0.460 0.551 + 0.0315i 0.281 

A-iii 11.34 0.460 1.15 

B 13.00 0.170 0.327 + 0.0126i 0.507 



fe f2 s : Surface angular velocity 
V: Eigenfrequency 

J ro C rt: Corotation radius; ro s : Surface radius 

mode is faster than the rotation of cylinder in all ra- 
dius, while the negative (A-i) corresponds to the op- 
posite. We also check the low T/|VT| bar instability 
in cylindrical model and find that the same behavior 
as m = 1 mode (B) appears in the distribution of the 
canonical angular momentum density (Fig. |SJ). 

We furthermore investigate m = 2 instability of 
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Figure 4: Distribution of canonical angular momentum density for m = 1 unstable mode (see Table ITT)l . Vertical 
dashed line represents the location of corotation radius of the mode. The Roman character in each panel corresponds 
to the model of the cylindrical gaseous fluid, respectively. Note that we normalized the distribution of the canonical 
angular momentum in an appropriate value, since the eigenfunction can be scaled arbitrarily. 
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Figure 5: Distribution of canonical angular momentum 
density for m — 2 unstable mode (see Table |HJ. The 
Roman character in the panel corresponds to the model 
of the cylindrical gaseous fluid. Vertical dashed lines 
mark the locations of corotation radius of the mode. 

Bardeen disk and classical bar instability of Maclaurin 
spheroid. In these cases, the canonical angular mo- 
mentum density that is analytically obtained is zero 
in all cylindrical radius. Therefore the behavior of the 
canonical angular momentum density shows a clear 
contrast between m = 1, 2 instability with high degree 
of differential rotation and m — 2 classical instability 
of uniformly rotating fluid. 

Finally, we adopt the method of canonical angular 
momentum to the nonlinear hydrodynamics. We iden- 
tify the complex eigenfrequency and the corotation 
radius from dipole or quadrupole diagnostics which is 
summarized in Table IIIII Note that we read off the 
eigenfrequency from those plots at the early stage of 
the evolution. The Eulerian perturbed velocity is de- 
fined by subtracting the velocity at equilibrium from 
the velocity. The Lagrangian displacement vector is 



Table III Eigenfrequency and the corotation radius of 3 
differentially rotating stars 



Model 


a [fi c ] 


ttcrt [Rcq] 


I 


0.590 + 0.0896i 


0.167 


II 


0.284 + 0.0121i 


0.492 


III 


0.757 + 0.200i 





extracted by using a linear formula for a dominant 
mode in each case. 

We show the snapshots of canonical angular mo- 
mentum density in Fig. El Since we determine the 
corotation radius using the extracted eigenfrequency 
and the angular velocity profile at equilibrium, the ra- 
dius does not change throughout the evolution. For 
low T/lVFl dynamical instability, the distribution of 
the canonical angular momentum drastically changes 
its sign around the corotation radius, and the maxi- 
mum amount of canonical angular momentum density 
increases at the early stage of evolution. This means 
that the angular momentum flows inside the corota- 
tion radius in the evolution. On the other hand, for 
high T/|VF| dynamical instability that is related to 
the classical bar instability, the distribution of the 
canonical angular momentum is smooth with no par- 
ticular feature and tends to have a positive portion 
outside. This means that the canonical angular mo- 
mentum flows outwards in the evolution, which is in 
clear contrast to the case of low T/IT4 7 ! one. 

From these different behaviors of the distribution 
of the canonical angular momentum, we see that the 
mechanisms working in the low T/|W| instabilities 
and the classical bar instability may be quite different, 
i.e., in the former the corotation resonance of modes 
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are essential, while the instability is global in the lat- 
ter case. 



4. Discussion 

We have studied the nature of dynamical one-armed 
spiral instability in differentially rotating stars both in 
linear eigenmode analysis and in hydrodynamic sim- 
ulation using canonical angular momentum distribu- 
tion. 

We have found that the one-armed spiral instabil- 
ity occurs around the corotation radius of the star by 
investigating the distribution of the canonical angu- 
lar momentum. We have also found by investigating 
the canonical angular momentum that the instability 
grows through the inflow of the angular momentum 
inside the corotation radius. The feature also holds 
for the dynamical bar instability in low T/|W|, which 
is in clear contrast to that of classical dynamical bar 
instability in high T/|W|. Therefore the existence of 
corotation point inside the star plays a significant role 
of exciting one-armed spiral mode and bar mode dy- 
namically in low T/|VU|. 

Finally, we mention the feature of gravitational 
waves generated from this instability. Quasi-periodic 
gravitational waves emitted by stars with m = 1 in- 
stabilities have smaller amplitudes than those emit- 
ted by stars unstable to the m = 2 bar mode. For 
m = 1 modes, the gravitational radiation is emitted 
not by the primary mode itself, but by the m = 2 sec- 
ondary harmonic which is simultaneously excited, pos- 
sibly through non-linear self-coupling of m=l mode. 
(Remarkably the precedent studies [jj, ||| found that 
the pattern speed of m = 2 mode is almost the same 
as that of m = 1 mode, which suggest the former is 
the harmonic of the latter.) Unlike the case for bar- 
unstable stars, the gravitational wave signal does not 
persist for many periods, but instead is damped fairly 
rapidly. 
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